Computation of soliton structure and analysis of chaotic behaviour in quantum deformed Sinh-Gordon model

Soliton dynamics and nonlinear phenomena in quantum deformation has been investigated through conformal time differential generalized form of q deformed Sinh-Gordon equation. The underlying equation has recently undergone substantial amount of research. In Phase 1, we employed modified auxiliary and new direct extended algebraic methods. Trigonometric, hyperbolic, exponential and rational solutions are successfully extracted using these techniques, coupled with the best possible constraint requirements implemented on parameters to ensure the existence of solutions. The findings, then, are represented by 2D, 3D and contour plots to highlight the various solitons’ propagation patterns such as kink-bright, bright, dark, bright-dark, kink, and kink-peakon solitons and solitary wave solutions. It is worth emphasizing that kink dark, dark peakon, dark and dark bright solitons have not been found earlier in literature. In phase 2, the underlying model is examined under various chaos detecting tools for example lyapunov exponents, multistability and time series analysis and bifurcation diagram. Chaotic behavior is investigated using various initial condition and novel results are obtained.


Introduction
Nonlinear models have earned considerable amount of attention due to their potential of demonstrating the dynamics of numerous natural and scientific phenomena, such as plasma waves (acoustic), gravitational waves, shallow water dynamics, fluid dynamics, nonlinear optics, and surface ocean waves.[1][2][3][4][5][6].Nonlinear partial differential equations (NLPDE) are essential for simulating exciting phenomenon in these areas and numerous real-world problems.Analytic solution of NLPDE, therefore, have been well researched, ever evolving and challenging field of exploration.Some eminent contributions are [7][8][9][10][11][12].
Quantum physics that regulates microscopic systems is a best example of the q deformed Sinh-Gordon equation [13].Numerous partial differential equations in classical calculus that have important applications in numerous fields have been studied by a number of scholars [14][15][16].These models play an instrumental role both in applied sciences and Mathematics.
With the advent of fractional calculus, it is now possible to explore mathematical models in a novel manner that incorporates fractional derivatives [17][18][19][20].Many scientists and researchers began integrating all the principles and relationships that are seen in fractional and traditional calculus into q calculus [21,22].Quantum calculus also referred as calculus without limits, is identical to conventional infinitesimal calculus but does not include the concept of limits.It provides ideas of h-calculus and q calculus where h refers Plank's constant and q representing quantum [23,24].Relativity theory, special functions, mathematical physics along with other areas have witnessed instrumental role of q calculus [25][26][27] and partial differential equations generated in q calculus are termed a q-deformed equations.
In this research, the following generalized q deformed Sinh-Gordon equation (Eleuch equation) [28] is taken into investigation.
Presently, theory of chaos and bifurcation is extensively used in examining differential equations.These are typically useful tools for understanding any intricate system that differential equations may control.A bifurcation is a qualitative shift driven by parameter change in the behaviour of a dynamical system.A variety of bifurcation schemes exist, including the saddle node bifurcation, Hopf bifurcations, period doubling bifurcation and pitchfork bifurcation.Examining the complex behaviour of nonlinear waves and investigating chaos theory are integral components of studying differential equation.The investigation of Chaos, a measure of stability when an external force is introduced to a nonlinear system, is essential in this contemporary era.Autonomous dynamic system's asymptotic behaviour is exclusively dictated by their introductory conditions.There are the four types of equilibrium behaviour: an equilibrium point, a limit cycle, a torus and chaos.
Recently, Jamal et al. [33] have studied the Novikov Veselov equation through bifurcation and chaos discovery tools and they also achieved the soliton solutions.Rafiq et al. [34] have investigated the shallow water waves through bifurcation and chaos analysis and also acquired the multi solitons.The conformable Fokas Lenells model has been examined by Lie and Huang [35], employing chaos theory and bifurcation analysis.Zhang et al. [36] carried out a investigation of bifurcation on the modified FitzHugh Nagumo neuron model and find out the novel results.Jamal et al. [37] have explored the nerve impulse model by practicing the phase portraits, quasi periodic, multistability, time series and sensitivity analysis, also obtained soliton solutions.
Liu and Li [38] examined the fractional perturbed Gerdjikov-Ivanov equation.Their research delved into the chaotic behavior of the model by introducing disturbance factors into the planar dynamical system.They analyzed various aspects of the model, including twodimensional and three-dimensional phase portraits, Poincare ´sections, and sensitivity analysis.Gu et al. [39] explored the (3+1)-dimensional negative-order Korteweg-de Vries-alogeroBogoyavlenskii-Schiff (KdV-CBS) equation, which extends the classical Korteweg-de Vries (KdV) equation and broadens the scope of nonlinear partial differential equations.
The aim of this article is to employ the new extended direct algebraic approach and modified auxiliary equation method to derive analytical solutions for the given model.As far as we are aware, this methods have not been used to evaluate Eq (1) in earlier research.The proposed approaches provides findings in several general and explicit form, covering trigonometric, exponential, rational and hyperbolic functions.This gives them various advantages over previously investigated techniques.Travelling wave transformation is used to convert nonlinear partial differential equation into nonlinear ordinary differential equation.The results are represented graphically as 3D and 2D charts for particular values of pertinent parameters.There are various methods for determining chaos [40].In this study, the most beneficial ones are emphasised.Here are Lyapunov exponents, multistability, time series and bifurcation diagram are discussed to review the chaotic behavior of the considered model.The study, in the opinion of the authors, is intriguing and has never been presented before the system in question.
The paper has been organised according to the format provided below.The mathematical framework and computation of soliton solution of the examined model have been described in Segment material and methods.Graphical illustrations and the chaotic behavior is investigated in Segment results and discussions.The research findings are highlighted in conclusion section.

The mathematical framework of the model
The subsequent transformation is proposed for the anatomization of the traveling wave solution of Eq (1): x À at ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi 1 À ðaÞ where the symbol α represents the speed of the travelling wave.Consequently, plugging Eq (3) into Eq (1), below mentioned ordinary differential equation has been derived.

Solutions by employing the extended direct algebraic method
In this segment, the generalized q deformed Sinh-Gordon equation is analyzed employing extended direct algebraic strategy.
where P is a polynomial involving higher order partial derivatives and nonlinear components in p(x, y, t).The travelling wave is transformed by a process known as: pðx; y; tÞ ¼ lnðmðZÞÞ; Z ¼ x þ y À at ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi 1 À ðaÞ 2 q ; ð8Þ here, α represents a nonzero constant.Eq (7) is transformed into a nonlinear ordinary differential equation of the following form: In Eq (9), the prime denotes the derivative with respect to η. Assume that Eq (9) has a solution of the following form: where while z, ν and μ are real constants, the value of Q can be determined by balancing the nonlinear terms and the highest order derivative appearing in Eq (9).
The value of Q is firstly obtained using homogeneous balance principle method.The value of Q = 2 is achieved by balancing the nonlinear terms with the highest order derivatives in Eq (6), which ultimately implies the solution in the known form: where P(η) satisfies Eq (11).Comparing the different powers of coefficients P(η), after plugging Eq (12) with Eq (11) into Eq (6) and have a system of suitable algebraic equations: The set of solutions below is produced by utilising a computational program to solve the aforementioned algebraic equations for the parameters a 0 , a 1 , a 2 , μ, ν, q and z: Case 1.If z 6 ¼ 0 and F < 0, then when the values of a 2 , μ and ν are entered via Eq (13) into Eq (12), the following results for Eq (6) are produced: Case 3. If ν = 0 and μz > 0, then Case 4. If ν = 0 and μz < 0, then Case 5.If μ = z and ν = 0, then Case 6.If z = −μ and ν = 0, then Case 7. If ν 6 ¼ 0 and μ = 0, then U 31 ¼ 4 ln ðrÞ 2 z 2 À mn zðcosh r ðnZÞ À sinh r ðnZÞ þ mÞ " # 2 ; Case 8.If ν = p, z = pq, (μ = 0 and q 6 ¼ 0), then

Solutions by employing the modified auxiliary equation method
In this section, the generalized q deformed Sinh-Gordon equation will be solved by using the modified auxiliary equation method.Consider Eq (9), which takes the following as its solution: where g(η) fits the auxiliary equation below and a i and b i are constants that must be computed.While where σ, β, θ and k are undefined constants with k > 0, k 6 ¼ 1.By balancing the nonlinear highest terms with highest derivatives into Eq (6), the value of n = 2 is obtained.The Eq (14) is convert in the following form as by using the above information: After employing Eq (15) with Eq (16) into Eq (6), set of algebraic equations turns up as follows and setting all the factors of the distinct powers of k g equal to zero, we have The set of solutions below is produced by utilising a computational program to solve the aforementioned algebraic equations for the parameters a 0 , a 1 , a 2 , b 1 , b 2 , β and σ: Plugging the values a 0 , a 1 , a 2 , b 1 ,b 2 , σ, and β via Eq (17) into Eq (16), it will turns up the examining solutions for the Eq (6): For β 2 − 4σθ < 0 and σ 6 ¼ 0 shows, ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi 4sy À b 2 p cot ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi 4sy À b 2 p cot ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi For β 2 − 4σθ > 0 and σ 6 ¼ 0 shows, ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi b 2 À 4sy p tanh ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi b 2 À 4sy p tanh ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi b 2 À 4sy p coth ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi b 2 À 4sy p coth ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi For β 2 − 4σθ = 0 and σ 6 ¼ 0 shows,

Investigation of chaotic behavior in perturbed dynamical system
This segment delves at the analysis of Eq (6), which describes quasi-periodic and chaotic behaviour.Employing a Galilean transformation, the studied equation is transformed into a planar dynamical system.To examine the chaotic behaviour of the planar dynamical system, a perturbation phrase O cos(θη) is introduced.Consequently, the dynamical planar system and the perturbation phrase have the following structure: with W = θη, is an independent system.The aforementioned system creates a disturbance phrase which represent the force and its frequency using the parameters O and θ.Behaviour of a system may change and looks chaotic whenever it is impacted by outside forces.In our research, we discover this behaviour in a chaotically behaving system (23) which demonstrate unpredictable time dependent trajectories that diverge from predictable trends.For determining chaos, we applied Lyapunov exponents method, multistability and time series analysis and bifurcation diagram.In order to comprehend the perturbed dynamical system, we then displayed the behaviour of the aforementioned exponents over time.We can figure out the chaotic behaviour of perturbed dynamical system at K 1 = 3, K 2 = −3.5,K 3 = 5.2, O = 3.7, θ = 1.8 and the starting conditions (0.6,0.6,0.6), (0.4,0.4,0.4) and (0.1,0.1,0.1)correspondingly by displaying the resulting Lyapunov exponents versus change in time in Figs (7)(8)(9).In (Fig 7), we have observed that the system ( 23) is chaotic at λ 1 = 0.0113, which is the largest positive Lyapunov exponent.This clearly shows that the orbits are diverging.In the same way, the largest positive Lyapunov exponent, λ 1 = 0.0313 for the perturbed dynamical system ( 23) is represents the presence of chaos.Moreover, the largest positive Lyapunov exponent, λ 1 = 0.0513 for the perturbed dynamical system ( 23) is represents the presence of chaos.
A characteristic of perturbed dynamical systems known as multistability and time series analysis denotes the presence of multiple alternative dynamic behaviors with an identical parameters set but diverse primary conditions.Among these behaviors are chaos, multistability, time series analysis periodicity and quasi periodicity which can manifest in the system in various situations.In Figs ( 10)-( 12), we have explored the multistability and time series examination of the perturbed system (23) under various beginning constraints.System (23) appears to be especially vulnerable to chaotic initial conditions based on observations.Comprehending the characteristics of multistability and time series analysis, which is an essential constituent of intricate dynamic systems, helps facilitate the explanation and prediction of these system's behaviours under diverse conditions.A visual representation that shows how a dynamical system behaves when a parameter is changed is known as bifurcation diagram.The magnitude or strength of the perturbation is the parameter of interest in the scenario of a dynamical system that is perturbed.One can learn more about the possible behaviours of the system, such as fixed points, chaos, or limit cycles, by looking at the bifurcation diagram.Bifurcation diagrams, especially, can be used to determine the key parameter values, such as the beginning of chaos or the shift from stable to unstable dynamics, at which the system experiences a qualitative change in behaviour.In (Fig

Conclusion
The q deformed Sinh-Gordon equation has been investigated by soliton structure and chaotic behaviour.The underlying equation is converted into an ordinary differential equation using the wave transformation technique.The new extended direct algebraic and modified auxiliary equation techniques have been effectively applied to obtain novel analytical traveling wave solutions.It is crucial to note that Sinh-Gordon equation solutions are achieved by utilizing hyperbolic, trigonometric, exponential, and rational functions.A variety of soliton structures have been produced for the resulting ordinary differential equation, taking into account different parameter values.Visual representations in both 2D and 3D forms are subsequently produced. ( The wave solution has been used to construct the many types of graphs.The structure of the travelling wave varies along with the unexplained variables of the outcomes.The composition of the solution has been examined.Now, the graphs are shown illustrating how the generalized q deformed Sinh-Gordon equation has changed over time.There are several different solution of the generalized q deformed Sinh-Gordon equation that contain unspecified parameters.These unknown factors have been effect on the nature of the solutions.In other words, if different specific values are assigned to the variables, several types of solutions can be produced from the general solution.The following diagram illustrates how the solution U 6 is affected by the parameters:The variables ρ, z, μ, ν, y, and α are all included in the inclusive finding U 6 .For the following values: ρ = 2, z = 1.5, μ = −0.5, ν = −1, y = 2.5, and α = 0.2 from solution U 6 , we acquire the kink-bright soliton structure.A 2D graph for oscillating with temporal components t = 0, 1 and 2 and components of velocity α = 0.2, 0.5, 0.8 within the limit −10 � x � 10 is shown in (Fig 1).Additionally, at intervals −5 � x � 5 and −5 � t � 5 and 3D graphs are displayed.Additionally, at intervals −5 � x � 5 and −5 � t � 5 and 3D graphs are displayed.The specific solution U 8 is made up of the parameters ρ, z, μ, ν, y and α.For the following values: ρ = 2, z = −1.8,μ = −0.5, ν = −2.4,y = −2 and α = 0.2 from solution U 8 , we are able to obtain the bright soliton structure.(Fig 2) shows a 2D graph for oscillating with temporal